Decoherence in Two Bose-Einstein Condensates 
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In this paper, decoherence in a system consisting of two Bose-Einstein condensates is investigated 
analytically. It is indicated that decoherence can be controlled through manipulating the interaction 
between the system and environment. The influence of the decoherence on quantum coherent atomic 
tunneling (AT) between two condensates with arbitrary initial states is studied in detail. Analytic 
expressions of the population difference (PD) and the AT current between two condensates are 
found. It is shown that the decoherence leads to the decay of the PD and the suppression of the AT 
current. 
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I. INTRODUCTION 



Recently, much attention has been paid to experimen- 
tal investigations [1-7] and theoretical studies [8-15] for 
systems consisting of two and multi Bose condensates 
since such systems give rise to a fascinating possibility 
of observing a rich set of new macroscopic quantum phe- 
nomena [16-20] which do not exist in a single condensate. 
Among important macroscopic quantum effects is the 
quantum coherent atomic tunneling (AT) between two 
trapped Bose condensates [15-18]. Several authors [19,20] 
showed that the AT can support macroscopic quantum 
self-trapping (MQST) due to the nonlinearity of atom- 
atom interactions in condensates. As is well known, no 
system can be completely isolated from its environment. 
In fact, in current experiments on trapped Bose con- 
densates of dilute alkali atomic gases condensated atoms 
continuously interact with non-condensate atoms (envi- 
ronment). Interactions between a quantum system and 
environment cause two types of irreversible effects: dis- 
sipation and decoherence [21,22]. Mathematically, the 
dissipation and decoherence can be understood in the 
following way. Let Hs and Hb be Hamiltonian of the 
system and environment (bath), respectively, and Hj be 
interacting Hamiltonian between the system and envi- 
ronment. When [Ht,Hs] ^ 0, which implies that the 
energy of the system is not conservative, the interaction 
Hi describes the dissipation. When [Hi, Hs] = 0, the en- 
ergy of the system is conservative, so the interaction Hi 
describes the decoherence. The dissipation effect, which 
dissipates the energy of the quantum system into the en- 
vironment, is characterized by the relaxation time scale 
T r . In contrast, the decoherence effect, which can be re- 
garded as a mechanism for enforcing classical behaviors 
in the macroscopic realm, is much more insidious because 
the coherence information leaks out into the environment 
in another time scale Td, which is much shorter than r r , 
as the quantum system evolves with time. Since macro- 
scopic quantum phenomena in Bose condensates mainly 



depend on rather than r r , the discussions in present 
paper only focus on the decoherence problem rather than 
the dissipation effect. To study the decoherence in Bose 
condensates is not only of theoretical interest but also 
importance from a practical point of view, since deco- 
hering would always be present in any Bose condensate 
experiments of trapped atoms. 

On the aspect of modeling dissipation and decoher- 
ence in trapped Bose condensates, some progress [9,23- 
26] has been made. In particular, Anglin [23] derived a 
master equation for a trapped Bose condensate by con- 
sidering a special model of a condensate confined in a 
deep but narrow spherical square-well potential. In his 
model the reservoir of non-condensate atoms consists of 
a continuum of unbound modes obtained by the scat- 
tering solutions of the potential well. Making use of the 
Anglin 's master equation, Ruostekoski and Walls [24] nu- 
merically simulated dissipative dynamics of a Bose con- 
densate in a double- well potential when the condensate is 
in the atomic coherent states, and shown that the inter- 
actions between condensate and non-condensate atoms 
make the MQST decay away. For a system consisting of 
two trapped weakly connected Bose condensates, there is 
quantum coherent AT between two condensates. Ques- 
tions that naturally arise are, what is the effect of deco- 
herence on the quantum coherent AT ? Does decoherence 
increase or decrease the AT current between them? In 
this paper, we analytically study the decoherence prob- 
lem in two Bose condensate, and investigate the influence 
of the decoherence on the quantum coherent AT between 
two trapped Bose condensates in terms of an exactly solv- 
able Hamiltonian. We will present analytic expressions 
of the population difference (PD) and the AT current 
between two Bose condensates, and show that the deco- 
herence leads to the PD decay and the suppression of the 
AT current. 

This paper is organized as follows. In Sec. II, we present 
an approximate analytic solution of the system consisting 
of two Bose condensates with a tunneling coupling with- 
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out decoherence. In Sec. Ill, we introduce a decoherence 
model and apply it to the two-condensate system, we dis- 
cuss the influence of decoherence the atomic tunneling. 
Concluding remarks are provided in the last section. 

II. TWO BOSE-EINSTEIN CONDENSATES 
WITH TUNNELING COUPLING 

Let us consider a system of two Bose condensates 
with weak nonlinear interatomic interactions and the 
Josephson-like coupling. Such a condensate system, in 
principle, can be produced in a double trap with two con- 
densates coupled by quantum tunneling and ground colli- 
sions, or in a system with two different magnetic sublevels 
of an atom, in which case the two species condensates cor- 
respond two electronic states involved. In the formalism 
of the second quantization, Hamiltonian of such a system 
can be written as 

H = H i + H 2 + Hi n t + Hj os , (1) 

H t = J dxV;i(x)[-|lv 2 + ^(x) 

+I7 i #(x)^(x)hfc(x),(* = l,2), (2) 

H int = U 12 f dx#(x)$(x)&(x)^(x), (3) 



H Jos =Aj dx[#(x)&(x) + &(x)$( x )]. (4) 

Here % — 1,2, ^j(x) and V>J(x) are the atomic field op- 
erators which annihilate and create atoms at position 
x, respectively. They satisfy the commutation relation 
[V'i(x), -0](x')] = SijS{x — x'). Hi and H 2 describe the 
evolution of each species in the absence of interspecies 
interactions. Hint describes interspecies collisions. Hj os 
is the Josephson-like tunneling coupling term. Atoms 
are confined in harmonic potentials V^(x)(i = 1,2). In- 
teractions between atoms are described by a nolinear 
self-interaction term Ui = 4-Kh 2 al c /m and a term that 
corresponds the nonlinear interaction between different 
species Ui 2 = ^irfi 2 a\ 2 /m, where a^ c is s-wave scattering 
lengths of species i and al 2 that between species 1 and 2. 
For simplicity, throughout this paper we set % = 1, and 
assume that a\ c = a 2 c = a sc , Vi(x) = V 2 (-x). 

It has been well known that the Hamiltonian (1) can 
reduce to a two-mode Hamiltonian [27] by the use of 
the approximation of the atomic field operators:^ (x) = 
Oi</>i(x), where a; = J dx0j(x)^j(x) are correspondent 
mode annihilation operators with real distribution func- 
tions (pi(x) and [a*, fit] = 1. Then the Hamiltonian (1) 
can be reduced to the two-mode Hamiltonian 



H = u)o(a\ai + a) 2 a 2 ) + q(a\ 2 a 2 + a^a 2 ) 

+g(a\a 2 + a\ai) + 2xa\aia\a 2 , (5) 

where q, \ and g are coupling constants which character- 
ize the strength of interatomic interaction in each con- 
densate, interspecies interaction, and the Josephson-like 
coupling, respectively. 

The valid conditions of the two-mode approxima- 
tion were demonstrated in Ref.[19] which indicate that 
the two-mode approximation is valid for weak many- 
body interactions, i.e., for small number of condensated 
atoms. As shown in Ref.[19], the two-mode approxi- 
mation should be acceptable for the number of atoms 
N < 2000 if the scattering length typically taken as 
a = 5nm for a large trap with the size L = lOfim. 

We note that the two-mode approximate Hamiltonian 
has the same form of that of a two-mode nonlinear opti- 
cal directional coupler [28]. It can not be exactly solved, 
but a closed analytical solution can be obtained under 
the rotating wave approximation suggested by Alodjanc 
et al. [29]. The approximate analytic solution is valid 
for the weak interactions between atoms, but it sheds 
considerable light on the AT under our consideration. 

In order to obtain an approximate analytic solution of 
the Hamiltonian (5), we introduce a new pair of bosonic 
operators: 

Ai = ±={ai + a 2 )e-^,A 2 = -^(oi - a 2 )e^ , (6) 

which satisfy the usual bosonic commutation relation: 
[Aj,Aj] = 5ij. Then the Hamiltonian (5) reduces to the 
following form 

H = QN + g(A\Ai - A\A 2 ) + ^q[(3N 2 - 2N) 

-{A\Ai-AlA 2 f] + l - x N 2 

- x A\AiA\A 2 + H', (7) 

where fl = (ujq — f ), the total number operator N — 
a\a,i + ala 2 = is a conservative constant, 

and H' is a nonresonant term which oscillates at the fre- 
quency Ag in the sense of the Alodjanc et aUs proposal 
[29]. The account of the fast oscillating terms results 
only in some additional oscillations which play no es- 
sential role in the evolution of the measurable quanti- 
ties specifying the macroscopic quantum phenomena of 
the two-condensate system, so that the nonresonant term 
are fully negligible. This is the rotating wave approxima- 
tion (RWA)in the sense of Ref.[29]. After neglecting the 
nonresonant term H', we get the following approximate 
Hamiltonian: 

H A = + g(A\Ai - A\A 2 ) + ^q[(3N 2 - 2N) 

-{A\Ai - A\A 2 f] + 1 - x n 2 - x a\AiA\A 2 . (8) 
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It is worthwhile noting that the dynamics of the non- 
RWA Hamiltonian (5) is often chaotic. A detailed investi- 
gation on chaotic behaviors in the two-condensate system 
is beyond the scope of present paper, and will be given 
elsewhere. Nevertheless, the RWA Hamiltonian (8) is an 
integrable Hamiltonian whose dynamics is regular, does 
not exhibit chaos. Hence, the the terms neglected in the 
RWA lead to chaos when they kept in the two-condensate 
system. This is very analogous to the case of the Jaynes- 
Cummings Model (JCM) [30] which describes the inter- 
action of a two-level atom with a single-mode electromag- 
netic field in quantum optics. It was well known that the 
RWA JCM is exactly solvable, but the non-RWA JCM 
[31] exhibits chaos. As shown in Ref.[31], the chaos was 
a consequence of inclusion of terms normally neglected 
in the RWA. 

Obviously, the Hamiltonian Ha is diagonal in the Fock 
space of the (Ai, A 2 ) representation defined by 



\n, to) 



1 



A\ n Al m \o,o), 



\Jn\m\ 



(9) 



where n and to take nonnegative integers. And we have 
HA\n, m) = E(n, m)\n, to) with the eigenvalues 

E{n, m) = (fi - |)(n + to) + g(n - to) + ^(3q + 2 X ) 
x(n + to) 2 - |(n - to) 2 - xnm. (10) 

For simplicity, all calculations below shall be car- 
ried out in the {A\, A 2 ) representation with the ba- 
sis: {|n, m),n,m = 0,1,2,...}, which is related to the 
(Si, 0,2) representation with a set of basis: {|?t.,to) = 
a[ n al m /Vn\m\\0, 0), n, m = 0,1,2,...} through the fol- 
lowing relation 



[n\m\{n — r + s)\(m — s + r)! 
2(«+™)/2( n - r )!( TO _ s )! 
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III. 



r=0 s=Q 

xe -i(2r-3m)7T/2| n _ r + Sj m _ s + r ). (11) 



INFLUENCE OF DECOHERENCE ON 
ATOMIC TUNNELING 



We now consider the effect of the decoherence. We use 
a reservoir consisting of an infinite set of harmonic os- 
cillators to model environment of condensate atoms in a 
trap, and we assume the total Hamiltonian to be 

H T = H A + J2^kb{b k + F({S})Y^ ck {b\ + k) 

k k 



where the second term is the Hamiltonian of the reser- 
voir. The last term in Eq.(12) is a renormalization term 



[32]. The third term in Eq.(12) represents the interac- 
tion between the system and the reservior with a cou- 
pling constant Ck, where {S} is a set of linear operators 
of the system or their linear combinations in the same 
picture as that of Ha, F({S}) is an operator function of 
{§}. In order to enable what the interaction between the 
system and environment describes is decoherence not dis- 
sipation, we require that the linear operator S commutes 
with the the Hamiltonian of the system H A - Then, the 
interaction term in Eq.(12) commutes with the Hamilto- 
nian of the system. This implies that there is no energy 
transfer between the system and its environment. So that 
it does describe the decoherence. The concrete form of 
the function FdS*}), which may be considered as an ex- 
perimentally determined quantity, may be different for 
different environment. Therefore, the decohering inter- 
action in (12) can not only describe decoherence caused 
by the effect of elastic collisions between condensate and 
non-condensate atoms for a Bosc condensate system, but 
also simulate decoherence caused by other decoherencing 
sources through properly choosing the operator function 
of the system F ({,§}). 

The Hamiltonian (12) can be exactly solved by making 
use of the following unitary transformation 



u = cx P [H A y2—(bt-bk)}. 



(13) 



Corresponding to the Hamiltonian (12), the total den- 
sity operator of the system plus reservoir can be ex- 
pressed as 



xe itJ2 k "kbl b *U e iHAt^ 



(14) 



In the derivation of the above solution, we have used 
p T (t) = U^p'^tp, where p' T = e - i6 ^ ^ T (0)e i6 ^ with 
H' T = UHtU- 1 and p' T {0) = Upr^Lf- 1 , where p T (0) 
the initial total density operator. 

We assume that the system and reservoir arc ini- 
tially in thermal equilibrium and uncorrelated, so that 
Pt(0) = p(0) ® pa, where p(0) is the initial density op- 
erator of the system, and pr the density operator of the 
reservoir, which can be written as pR = Y\ k pk{^) with 
/ifc(O) is the density operator of the fc-th harmonic oscil- 
lator in thermal equilibrium. After taking the trace over 
the reservoir, from Eq.(14) we can get the reduced den- 
sity operator of the system, denoted by p(t) = trRprit), 
its matrix elements in the {A\, A 2 ) representation arc 
explicitly written as 



(0)R(m',n')( m ,n) (t) 



xe 



-i[F({S(m' ,n')})-F({S(m,n)})]t 



(15) 
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where F({S(m, n)}) is an eigenvalue of the operator func- Q 2 {t) = 0, then we get a expression of the PD between 
tion FdS 1 }) in an eigenstate of Ha- R( m ',n')(m,n){t) is a two condensates 
reservoir-dependent quantity given by 

K*) = -2^^V / sF+l)|P(r+l, S -l)(r >S )(0)| 



R 



(m' in'){m,n) 



(*) 



xe 



t [_F 2 ({S(m'V)})--F 2 ({S(m,r l )})]Q 1 (t) 
-{F({S(m',n')})-F({S(m,n)})] 2 Q 2 (t) 



(16) 



where the two reservoir-dependent functions are given by 



Qi(t) 



dujJ((jj) 



c 2 n 



0J< 



sin(wt), 



Qi{t) 



<hjJ{ui) 



c 2 (c) 







sin (_)coth(— ), 



(17) 



(18) 



Here we have taken the continuum limit of the reservoir 
modes: J2 k — ► J °°dwJ(w), where J(w) is the spectral 
density of the reservoir, c(w) is the correspondeing con- 
tinuum expression for Cfe, and /3 = 1/fcgT with fc^ and T 
being the Boltzmann constant and temperature, respec- 
tively. 

It is well known that decoherence corresponds to the 
decay of off-diagonal elements of the reduced density ma- 
trix of a quantum system. For the case under our consid- 
eration, the degree of decoherence is determined by the 
decaying factor in Eq.(16). It is interesting to note that 
if we choose a proper operator function F({S}) to make 
F({S(m',n')}) = F({S(m,n)}) for (m',n') ^ (m,n), 
then we find that 



P{m' ,n'){m,n)\t) P{m' ,n')(m,n) (0) 



(19) 



which indicates that the quantum system maintains its 
initial quantum coherence, namely, the time evolution 
of decoherence-free of the quantum system is realized. 
Therefore, we conclude that one can control decoherence 
by manipulating interaction function i^dS 1 }). 

Eqs.(15) and (16) indicate that the interaction between 
the system and its environment induces a phase shift and 
a decaying factor in the reduced density operator of the 
system. We now consider the PD between the two con- 
densates in the presence of the decoherence, defined by 
p(t) = Nt(t) - N 2 (t) with Ni = (a\a,i). We find that 



P( f ) = - 2 E E + l)|P(r+l,.-l)(r,a)(0)| 

r s 

x sm[9 rs - v-(t - v+ s Qi(t))]e~ {v ™ fQ2{t) , (20) 



where we have introduced the symbols: 

P(r+l,s-l)(r,s)(0) = \P(r+l,s-l)(r,s)(0)\e ters , 

w± = F({S(r + l,a- 1)}) ± F({S(r, *)}). 



(21) 



(22) 



From Eq.(20) we see that if we do not take into ac- 
count the influence of the decoherence, i.e., set Qi{t) = 



x sin(# rs - v rs t), 



(23) 



which implies that the time evolution of the PD is pe- 
riodic. In particular, if we take ^({5*} = Ha, we find 
that when 2g/(q — x) = K { being an integer), we have 
a nonzero time-average value of the PD 



P = -2 ^ V(r-K+ l)(r + l)\p {r +Ur-K)^r-K + l)(0)\ 

r 

x sin(0 rr _ K+ i), (24) 

which means that the two-condensate system under our 
consideration exhibits the MQST [18,19] when the deco- 
herence is absent. 

The coherent AT current between the two condensates, 
defined by I(t) = N^t) - N 2 (t), is given by 

I(t) = 2 ^ ^ V*(r + l)«-|p( r+ i, a _i)( r ,a) (0)| 

r s 

x{(l - v+Qi{t)) cos[6 rs - v~(t - v+Mt))] 
+v~ s Q2(t) sin[6» rs - v~ s (t - v+ s Qi(i))]} 

xe -(v- s ) 2 Q2(t)_ (25) 

From Eqs.(18), (20) and (25) we can immediately draw 
one important qualitative conclusion: since Q 2 (t) is pos- 
itive definite, the existence of the decoherence is always 
to tend to suppress the PD and the AT current between 
the two condensates. This answers the question: "Does 
the decoherence increase or decrease the AT?" . 

From Eqs.(17),(18), (20) and (25) we sec that all nec- 
essary information about the effects of the environment 
on the PD and the AT current is contained in the spec- 
tral density of the reservoir. To proceed further let us 
now specialize to the Ohmic case [33] with the spectral 
distribution J(uj) = t i(^ e~"/"'=, where lo c is the high 
frequency cut-off, r\ is a positive characteristic parameter 
of the reservoir. With this choice, at low temperature 
the functions Qi(t) and Q 2 {t) are given by the following 
expressions 



Qi(t) = 77 tan 1 (u> c t), 



(26) 



Q 2 (t) = r,{\ ln[l + M) 2 ] + ln[A sinh(^)]}. (27) 

Recent experiments [1,7] on two condensates have es- 
tablished a typical time scale at which the two conden- 
sates preserve coherence. The value of the typical time 
scale is t^lOOms. In the meaningful domain of time 
uj c t S> 1 which requires uj c S> 10Hz which can be easily 



4 



satisfied for an usual reservoir [34,32,33], at zero temper- 
ature, we have Qi(t) = r]/(uj c t 2 ), and Q2(t) = i]ln(ijj c t), 
then we find 



P(t) = \AF+l)|p(r+l,s-l)(r, S )(0)| 
r s 

x sm[9 rs ~ v~(jt - v+ s Qx{t))]^ c t)-^ v rs)\ ( 28 ) 

X {(1 - 5H«t- a ) cos[0 rs - «- (* - «+Qi(t))] 

+»7tV s *~ 1 sin[0 r;5 - u~(f - u^Qi (*))]} 
x(w c t)-"(««) a , (29) 

which indicate that the PD and the AT current decay 
away according to the "power law" , where we have noted 
that the decaying factors can not be taken outside the 
summation on the r.h.s. of Eqs.(28) and (29). 

At finite temperature, we have Qi(t) = i]/(uj c t 2 ), and 

Q2(t) = v[H^r) + f], so that 



*>(*) = - 2 E E + l)IP(r+l,.-l)(r, S )(0)| 

r s 

x sin[6» rs - v~ s (t - w+Qx(t))] 
x(^)-"^) 2 ex P [-^^i], (30) 

J W = 2 EE^ + l)^|P(r+l, a -l)(r,«)(0)| 



the PD and the AT current between two condensates. 
We have obtained analytic expressions of the PD and 
the AT current and found that for the reservoir-spectral 
density of the Ohmic case, the PD and the AT current 
decay away by the "power law" at zero temperature; at 
finite temperature, the PD decays away by the "expo- 
nential law" while the decay of the AT current contains 
both "exponential-law" and "power-law" components. It 
is worthwhile to note that our results are obtained for 
arbitrary initial states of the two condensates, our entire 
analysis is carried out without invoking the assumption of 
Bose-broken symmetry which has recently been shown to 
be unnecessary for a Bose condensate of trapped atoms 
[35]. Also it should be pointed out that these results 
are obtained under the RWA in the sense of Alodjanic et 
aVs proposal, they are valid for interatomic weak non- 
linear interactions in two condensates. The RWA essen- 
tially changes the two-condensate system into an inte- 
grable system, hence it suppresses the chaotic behaviors 
of the two-condensate system. 
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X {(i - a^ t - 2 ) cos[e rs - v-(t - V+Mt))] 
+T^s. sin[9rs _ „- { t v+Mt))]} 

x(^)-(0 2 ex P [-I^i], (31) 

which indicate that at finite temperature the PD decays 
away according to the "exponential law" , and the decay 
of the AT current becomes more complicated than that of 
the PD due to the factor {j]V^ s /uj c )t~ 2 before the cosine 
function in Eq.(31). 



IV. CONCLUDING REMARKS 

We have present a decoherence model which is exactly 
solvable, and applied it to study decoherence in two Bose 
condensates. We have indicated that one can control de- 
coherence by manipulating interaction between a quan- 
tum system and environment. We have investigated the 
influence of the decoherence on quantum coherent AT be- 
tween two trapped Bose condensates with arbitrary ini- 
tial states, and shown that the decoherence suppresses 
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